A&A  479,  L21-L24  (2008) 

DOI:  10.1051/0004-6361:20079237 
©  ESO  2008 


Astronomy 

Astrophysics 


Letter  to  the  Editor 


J,  H,  K  spectro-interferometry  of  the  Mira  variable  S  Orionis* 

M.  Wittkowski1,  D.  A.  Boboltz2,  T.  Driebe3,  J.-B.  Le  Bouquin4  F.  Millour3  K.  Ohnaka3,  and  M.  Scholz56 

1  ESO,  Karl-Schwarzschild-Str.  2,  85748  Garching  bei  Miinchen,  Germany 
e-mail:  mwittkow@eso .  org 

2  US  Naval  Observatory,  3450  Massachusetts  Avenue,  NW,  Washington,  DC  20392-5420,  USA 

3  Max-Planck-Institut  fur  Radioastronomie,  Auf  dem  Hiigel  69,  53121  Bonn,  Germany 

4  ESO,  Casilla  19001,  Santiago  19,  Chile 

5  Institut  fur  Theoretische  Astrophysik  der  Univ.  Heidelberg,  Albert-Ueberle-Str.  2,  69120  Heidelberg,  Germany 

6  Institute  of  Astronomy,  School  of  Physics,  University  of  Sydney,  Sydney  NSW  2006,  Australia 

Received  12  December  2007/  Accepted  29  December  2007 

ABSTRACT 

Aims.  We  present  J,  H,  K  spectrally  dispersed  interferometry  with  a  spectral  resolution  of  35  for  the  Mira  variable  S  Orionis.  We  aim 
at  measuring  the  diameter  variation  as  a  function  of  wavelength  that  is  expected  due  to  molecular  layers  lying  above  the  continuum¬ 
forming  photosphere.  Our  final  goal  is  a  better  understanding  of  the  pulsating  atmosphere  and  its  role  in  the  mass-loss  process. 
Methods.  Visibility  data  of  S  Ori  were  obtained  at  phase  0.78  with  the  VLTI/ AMBER  instrument  using  the  fringe  tracker  FINITO 
at  29  spectral  channels  between  1.29yum  and  2.32/rm.  Apparent  uniform  disk  (UD)  diameters  were  computed  for  each  spectral 
channel.  In  addition,  the  visibility  data  were  directly  compared  to  predictions  by  recent  self-excited  dynamic  model  atmospheres. 
Results.  S  Ori  shows  significant  variations  in  the  visibility  values  as  a  function  of  spectral  channel  that  can  only  be  described  by  a  clear 
variation  in  the  apparent  angular  size  with  wavelength.  The  closure  phase  values  are  close  to  zero  at  all  spectral  channels,  indicating 
the  absence  of  asymmetric  intensity  features.  The  apparent  UD  angular  diameter  is  smallest  at  about  1.3  /im  and  1.7  fim  and  increases 
by  a  factor  of  ~1.4  around  2.0/rm.  The  minimum  UD  angular  diameter  at  near-continuum  wavelengths  is  ©UD  =  8.1  ±  0.5  mas, 
corresponding  to  R  ~  420  RB.  The  S  Ori  visibility  data  and  the  apparent  UD  variations  can  be  explained  reasonably  well  by  a  dynamic 
atmosphere  model  that  includes  molecular  layers,  particularly  water  vapor  and  CO.  The  best-fitting  photospheric  angular  diameter  of 
the  model  atmosphere  is  ©phot  =  8.3  ±  0.2  mas,  consistent  with  the  UD  diameter  measured  at  near-continuum  wavelengths. 
Conclusions.  The  measured  visibility  and  UD  diameter  variations  with  wavelength  resemble  and  generally  confirm  the  predictions 
by  recent  dynamic  model  atmospheres.  These  size  variations  with  wavelength  can  be  understood  as  the  effects  from  water  vapor 
and  CO  layers  lying  above  the  continuum-forming  photosphere.  The  major  remaining  differences  between  observations  and  model 
prediction  are  very  likely  due  to  an  imperfect  match  of  the  phase  and  cycle  combination  between  observation  and  available  models. 

Key  words,  techniques:  interferometric  -  stars:  AGB  and  post-AGB  -  stars:  atmospheres  -  stars:  individual:  S  Ori  -  stars:  mass-loss 


1.  Introduction 

Mira  stars  are  low-mass,  large-amplitude,  long-period  variable 
stars  on  the  AGB,  evolving  toward  the  planetary  nebula  and 
white  dwarf  phases.  They  exhibit  a  mass-loss  rate  on  the  order 
of  ~  1CT6  Mo/year  that  significantly  affects  the  further  stellar  evo¬ 
lution  and  is  one  of  the  most  important  sources  for  the  chemical 
enrichment  of  the  interstellar  medium.  The  dust  condensation  se¬ 
quence,  the  wind-driving  mechanism,  and  the  role  of  pulsation 
are  currently  not  well  understood,  in  particular  for  oxygen-rich 
AGB  stars  (Woitke  et  al.  2006;  Hofner  &  Andersen  2007).  The 
pulsating  atmospheres  of  Mira  stars  can  become  very  extended 
because  of  dynamic  effects  including  shock  fronts,  and  they  are 
very  cool  in  their  outer  parts.  Here,  molecules  can  form,  which 
for  O-rich  stars  are  most  importantly  H2O,  CO,  TiO,  and  SiO 
(Tsuji  et  al.  1997;  Tej  et  al.  2003;  Ohnaka  2004).  Wittkowski 
et  al.  (2007)  found  for  the  case  of  S  Ori  that  AI2O3  dust 
condenses  within  the  extended  atmosphere  at  phase-dependent 
distances  of  1. 8-2.4  photospheric  radii.  This  extended  atmo¬ 
sphere,  which  is  characterized  by  phase-dependent  temperature 
and  density  stratifications,  the  presence  of  molecular  layers, 
and  the  formation  of  dust,  is  thus  of  particular  interest  for  our 

*  Based  on  observations  made  with  the  VLT  Interferometer  (VLTI)  at 
Paranal  Observatory  under  program  ID  080.D-0691. 


better  understanding  of  pulsation  and  mass  loss.  Observed  radii 
of  Mira  stars  have  been  found  to  differ  for  different  optical 
and  infrared  bandpasses  (e.g.  Thompson  et  al.  2002;  Mennesson 
et  al.  2002;  Ireland  et  al.  2004a;  Perrin  et  al.  2004;  Eisner  et  al. 
2007),  and  this  has  been  attributed  to  the  presence  of  molec¬ 
ular  layers  located  above  the  continuum-forming  photosphere. 
Here,  we  present  both  a  spectro-interferometric  observation  of 
the  Mira  star  S  Ori  that  covers  the  near-infrared  J,  H,  and  K 
bands  simultaneously  at  a  spectral  resolution  of  35  and  a  com¬ 
parison  to  recent  self-excited  dynamic  model  atmospheres. 

S  Ori  is  a  Mira  variable  star  with  spectral  type  M6.5e-M9.5e 
and  V  magnitude  7.2-14.0  (Samus  et  al.  2004).  We  use  a  Julian 
Date  of  last  maximum  brightness  Tq  —  2  453  190  days,  a  period 
P  =  430  days  (as  in  Wittkowski  et  al.  2007),  and  the  distance 
of  480  pc  ±  120  pc  from  van  Belle  et  al.  (2002).  The  broad¬ 
band  near-infrared  K  UD  angular  diameter  of  S  Ori  has  been 
measured  by  van  Belle  et  al.  (1996),  Millan-Gabet  et  al.  (2005), 
and  Boboltz  &  Wittkowski  (2005)  to  values  between  9.6  mas  and 
10.5  mas  at  different  phases.  Joint  VLTI/MIDI  and  VLBA/SiO 
maser  observations  by  Wittkowski  et  al.  (2007)  have  shown  that 
S  Ori  exhibits  significant  phase-dependencies  of  the  atmospheric 
extension  and  dust  shell  parameters  with  photospheric  angular 
diameters  between  7.9  mas  and  9.7  mas. 
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Table  1.  Observation  log.  Night  starting  12  October  2007,  JD  2454  386. 


Target 

Purpose 

Old 

[mas] 

DIT 

[ms] 

Time 

[UTC] 

®Vis 

Sp  [m] 

E0-G0/G0-H0/E0-H0 

PAp 

deg 

AM 

Seeing 

["] 

To 

[ms] 

45  Eri 

Calibrator  (K3  Il-III) 

2.15  ±0.04 

25 

08:03-08:07 

16.0/31.9/47.9 

-107 

1.1 

1.3 

1.3 

45  Eri 

Calibrator  (K3  II-III) 

2.15  ±0.04 

50 

08:09-08:12 

16.0/32.0/48.0 

-107 

1.1 

1.3 

1.3 

y  Eri 

Check  star  (M0.5  Illb) 

8.74  ±0.09 

25 

08:22-08:26 

15.6/31.2/46.8 

-104 

1.1 

1.4 

1.2 

y  Eri 

Check  star  (M0.5  Illb) 

8.74  ±  0.09 

50 

08:28-08:32 

15.5/31.0/46.5 

-104 

1.1 

1.4 

1.2 

SOri 

Science  target 

25 

08:45-08:49 

0.78 

15.9/31.8/47.8 

-107 

1.1 

1.3 

1.3 

SOri 

Science  target 

50 

08:52-08:57 

0.78 

16.0/31.9/47.9 

-107 

1.1 

1.3 

1.3 

a  Hor 

Calibrator  (K2  III) 

2.76  ±  0.03 

25 

09:14-09:18 

14.5/28.9/43.3 

-91 

1.1 

1.2 

1.4 

a  Hor 

Calibrator  (K2  III) 

2.76  ±0.03 

50 

09:21-09-24 

14.3/28.7/43.0 

-90 

1.1 

1.4 

1.2 

SOri 

Science  target 

25 

09:36-09:40 

0.78 

15.9/31.8/47.7 

-106 

1.1 

1.5 

1.1 

Wavelength  (nm)  Wavelength  (nm) 


Fig.  1.  Measured  S  Ori  visibility  data  compared  to  models  of  a  UD  with  a  constant  diameter  (red  dashed  lines),  of  a  Gaussian  disk  of  constant 
diameter  (green  dashed  line),  and  of  the  M18n  atmosphere  model  (blue  solid  line).  For  the  projected  baseline  lengths  and  angles  see  Table  1. 


Wavelength  (nm)  Wavelength  (nm)  Wavelength  (nm) 


Wavelength  (nm) 


Fig.  2.  Measured  y  Eri  visibility  data  compared  to  models  of  a  UD  with  a  constant  diameter  (red  dashed  line)  and  of  an  ATLAS  9  model  atmosphere 
with  Teff  =  3750  K,  logg  =  1.5,  solar  chemical  abundance  (blue  solid  line).  For  the  projected  baseline  lengths  and  angles  see  Table  1. 


2.  Observations  and  data  reduction 

We  obtained  near-infrared  J,  H,  K  interferometry  of  S  Ori  with 
the  instrument  AMBER  (Petrov  et  al.  2007)  in  low-resolution 
mode  at  the  ESO  VLTI  using  the  fringe  tracker  F1N1TO  and 
three  VLTI  Auxiliary  Telescopes  (ATs)  on  12  October  2007 
(JD  2454  386).  The  ATs  were  positioned  on  stations  E0,  GO, 
and  HO.  The  date  of  observation  corresponds  to  a  visual  phase 
Ovis  =  0.78,  with  an  uncertainty  of  about  0.1.  The  details  of 
the  observing  sequence  are  listed  in  Table  1,  including  the  pro¬ 
jected  baseline  lengths  ( Bp )  and  position  angles  (PAp),  the  air- 
mass  (AM),  and  the  optical  seeing  and  coherence  time.  Ambient 
conditions  were  not  very  good  but  stable.  The  airmass  was  the 
same  for  all  observations.  Data  were  recorded  using  two  differ¬ 
ent  detector  integration  times  (DITs)  of  25  ms  and  50  ms. 

The  VLTI  fringe-tracker  FINITO  records  fringes  on  the  two 
shortest  baselines  using  70%  of  the  H- band  light.  Output  signals 
are  processed  in  realtime  and  used  to  compensate  for  the  fringe 
motion  due  to  atmospheric  turbulence.  Owing  to  the  low  coher¬ 
ence  time  during  the  observations,  FINITO  was  only  able  to  pro¬ 
vide  average  performances  of  0.2-0. 4  /im  RMS  (to  be  compared 
with  0.1  /rm  RMS  achieved  in  good  conditions).  Nevertheless 
such  a  performance  is  sufficient  for  increasing  the  signal-to- 
noise  ratio  (S/R)  of  the  AMBER  data  and  stabilizing  the  transfer 
function.  All  observations  used  here  were  obtained  with  exactly 


the  same  FINITO  controller  parameters,  which  is  important  for 
avoiding  systematic  biases  in  calibrating  the  absolute  visibility. 

In  addition  to  S  Ori,  the  two  calibration  stars  45  Eri  and 
a  Hor  were  observed  close  in  time,  as  was  y  Eri.  The  last  is  a 
regular  non-pulsating  M  0.5  giant  with  a  well-known  angular  di¬ 
ameter  similar  to  that  of  S  Ori  which  is  not  expected  to  exhibit 
strong  effects  from  molecular  layers.  This  data  is  used  to  check 
that  any  strong  wavelength-dependent  features  found  for  S  Ori 
are  not  caused  by  any  systematic  effects  of  the  instrument.  The 
spectral  types  and  angular  diameters  of  the  calibration  and  check 
stars  are  from  Borde  et  al.  (2002). 

Raw  visibility  and  closure  phase  values  were  computed 
using  the  latest  version  of  the  amdlib  package  (version  2.0 
beta  2b)  and  the  yorick  interface  provided  by  the  AMBER 
consortium  and  the  Jean-Marie  Mariotti  Center.  Absolute  wave¬ 
length  calibration  was  performed  by  correlating  the  raw  spectra 
of  all  four  stars  with  a  model  of  the  atmospheric  transmission 
with  the  same  spectral  resolution.  In  particular,  we  used  a  plateau 
in  the  transmission  curve  at  A  ~  2.0 /nil.  The  offset  with  respect 
to  the  original  wavelength  table  was  3  spectral  channels  (0. 1  /nm 
at  a  wavelength  of  2.0 /nm).  We  estimated  the  error  of  the  ab¬ 
solute  wavelength  calibration  to  1  pixel  (-0.03  /nm).  Individual 
frames  were  averaged  after  frame  selection  keeping  70%  of  the 
best  frames  based  on  piston  (to  remove  the  frames  when  the 
FINITO  loop  was  not  closed)  and  out  of  these  keeping  30%  of 
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Wavelength  (nm) 

Fig.  3.  S  Ori  UD  diameter  values  as  a  function  of  wavelength  compared 
to  the  prediction  by  the  M18n  model  atmosphere.  Also  indicated  are  the 
positions  of  H2O  and  CO  bands  after  Lan§on  &  Wood  (2000). 

the  best  frames  based  on  S/R.  We  verified  that  keeping  up  to  80% 
of  the  best  frames  based  on  S/R  did  not  significantly  change  the 
results.  The  S  Ori  and  y  Eri  visibility  data  were  calibrated  for 
each  DIT  value  separately  using  the  45  Eri  and  a  Hor  calibra¬ 
tion  star  data.  After  calibration,  the  different  calibrated  S  Ori 
and  y  Eri  data  were  averaged.  The  errors  of  the  calibrated  vis¬ 
ibility  data  include  the  statistical  error  of  averaging  the  single 
frames,  the  errors  of  the  calibration  stars’  angular  diameters,  and 
the  variation  of  the  available  transfer  function  measurements. 


3.  Results 

Figures  1  and  2  show  the  resulting  visibility  and  closure  phase 
data  for  S  Ori  and  for  the  check  star  y  Eri,  respectively.  The  gaps 
in  the  visibility  data  around  1.45  gm  and  1.85  fjm  correspond  to 
the  regions  between  the  bands.  Also  shown  are  the  best  fitting 
models  of  a  UD  with  a  constant  diameter,  of  a  Gaussian  disk 
with  a  constant  diameter,  and  of  atmosphere  models.  The  lat¬ 
ter  are  for  S  Ori  the  M18n  model  (described  in  detail  below  in 
Sect.  4),  and  for  y  Eri  an  ATLAS  9  model  atmosphere  (Kurucz 
1993)  with  Teff  =  3750  K,  log  g  =  1.5,  and  solar  chemical  abun¬ 
dance  as  expected  for  y  Eri  (Borde  et  al.  2002).  The  comparison 
of  the  S  Ori  visibility  data  to  the  dynamic  model  atmosphere  is 
described  in  detail  below  in  Sect.  4. 

The  calculation  of  synthetic  visibility  values  and  the  fits  to 
the  interferometric  data  were  performed  as  in  Wittkowski  et  al. 
(2006,  2007).  The  fitted  angular  diameter  of  the  Gaussian  model 
corresponds  to  the  FWHM,  that  of  the  plane-parallel  ATLAS  9 
model  to  the  0%  intensity  (limb-darkened)  radius,  and  that  of  the 
M18n  model  to  the  1.04 /urn  (photospheric)  radius  (as  defined  in 
Ireland  et  al.  2004c).  The  best-fit  angular  diameters  are  S  Ori: 
0UD  =  10.8  mas;  ©Gaussian  =  6.9  mas;  ©Mi8n  =  8.3  mas;  y  Eri: 
©ud  =  8.5  mas;  ©atlas9  =  8.9  mas.  The  errors  are  <x  ~  0.2  mas. 

The  visibility  data  of  the  check  star  y  Eri  can  be  described 
well  by  a  UD  of  constant  diameter  and  by  the  ATLAS  9  model 
atmosphere.  There  are  no  significant  wavelength-dependent  de¬ 
viations  between  measured  visibility  data  and  the  UD  model.  It 
is  not  yet  clear  whether  the  relatively  low  ./-band  visibilities  for 
the  G0-H0  baseline  and  the  deviations  in  the  closure  phase  val¬ 
ues  near  the  flip  are  caused  by  an  asymmetric  stellar  surface  fea¬ 
ture  or  a  systematic  calibration  uncertainty.  Note  that  the  angular 
diameter  of  y  Eri  based  on  a  given  model  is  well-constrained  by 
the  position  of  the  visibility  minimum,  which  is  independent  of 
an  absolute  visibility  calibration.  The  resulting  angular  diameter 
©atlas9  =  8.9  ±  0.2  mas  is  consistent  with  the  value  given  in 
Borde  et  al.  (2002)  of  ©ld  =  8.74  ±  0.09  mas. 


The  visibility  data  of  S  Ori  show  significant  wavelength- 
dependent  features  clearly  deviating  from  UD  and  Gaussian 
models  of  constant  diameter  on  all  three  baselines.  This 
indicates  variations  in  the  apparent  angular  diameter.  The  S  Ori 
closure  phase  values  are  consistent  with  zero  at  all  spectral 
channels,  indicating  the  absence  of  asymmetric  features  in 
the  intensity  distribution.  In  order  to  characterize  the  variation 
of  S  Ori’s  angular  diameter  as  a  function  of  wavelength, 
we  fitted  UD  diameters  to  the  data  of  each  spectral  channel 
separately.  Figure  3  shows  the  resulting  UD  diameter  values 
as  a  function  of  wavelength.  Note  that  the  intensity  profile  of 
a  Mira  star  is  generally  not  expected  to  be  close  to  a  UD  and 
that  this  approach  can  only  give  a  rough  estimate  of  S  Ori’s 
characteristic  size  as  a  function  of  wavelength.  Fits  of  Gaussian 
functions  lead  to  similar  results.  The  apparent  UD  angular 
diameter  shows  clear  variations  with  wavelength.  It  is  smallest 
at  about  1 .3  //m  and  1 .7  /mi  and  increases  by  a  factor  of  ~1.4 
around  2.0 /jm.  The  minimum  UD  angular  diameter  of  S  Ori 
at  the  near-continuum  wavelengths  is  ©ud  =  8.1  ±  0.5  mas. 
This  result  at  visual  phase  0.78  is  consistent  with  the  S  Ori 
photospheric  angular  diameters  between  7.9  mas  at  phase 
0.55  and  9.7  mas  at  phase  1.15  derived  in  Wittkowski  et  al. 
(2007)  based  on  VLTI/MIDI  data  and  modeling  with  dynamic 
model  atmospheres  and  a  radiative  transfer  model  of  the 
dust  shell.  These  photospheric  angular  diameters  are  also 
consistent  with  previous  broadband  measurements  (cf. 
Wittkowski  et  al.  2007).  With  the  adopted  distance  to  S  Ori,  the 
angular  diameter  ©ud  =  8.1  ±  0.5  mas  corresponds  to  a  radius 
418  ±  130Ro. 

4.  Comparison  to  dynamic  Mira  star  atmosphere 
models 

Few  dynamic  atmosphere  models  for  oxygen-rich  Mira  stars  are 
available  that  include  the  effects  of  molecular  layers.  The  P  and 
M  model  series  (Ireland  et  al.  2004b, c)  are  complete  self-excited 
dynamic  model  atmospheres  of  Mira  stars  designed  to  match  the 
prototype  oxygen-rich  Mira  stars  o  Cet  and  R  Leo.  They  have 
been  used  successfully  for  comparisons  to  recent  broadband  in¬ 
terferometric  data  of  o  Cet  and  R  Leo  (Woodruff  et  al.  2004; 
Fedele  et  al.  2005).  Compared  to  o  Cet  and  R  Leo,  S  Ori  is  a 
slightly  cooler  Mira  variable  with  a  longer  period,  a  higher  main- 
sequence  precursor  mass,  and  a  larger  radius.  However,  when 
scaled  to  variability  phases  between  0  and  1  and  to  the  corre¬ 
sponding  angular  size  on  the  sky,  the  general  model  results  are 
not  expected  to  be  dramatically  different  for  S  Ori  compared  to 
o  Cet  and  R  Leo  (cf.  the  discussion  in  Wittkowski  et  al.  2007). 
The  M  model  series  was  chosen  to  model  the  atmosphere  of 
S  Ori  by  Wittkowski  et  al.  (2007)  as  the  currently  best  avail¬ 
able  option  to  describe  Mira  star  atmospheres.  Monochromatic 
center-to-limb  variations  (CLVs)  at  46  angles  between  0  and  5  Rp 
based  on  the  P  and  M  models  were  recomputed  for  the  wave¬ 
length  range  from  1-2.5  /mi  in  steps  of  0.001  /mi. 

The  P  and  M  dynamic  model  atmospheres  predict  signifi¬ 
cant  changes  in  the  monochromatic  radius  R,\  =  R,i(t,i  =  1), 
caused  by  molecular  layers  that  lie  above  the  continuum-forming 
photosphere  and  significantly  affect  certain  bandpasses.  Figure  4 
shows  the  monochromatic  radius  R,\  =  R,{(t/1  =  1)  in  units 
of  the  non-pulsating  parent  star  radius  for  the  example  of  three 
models  of  the  series.  It  illustrates  the  strong  phase  dependence 
of  the  molecular  layers.  The  models  are  M16n  (model  phase 
1.60),  M18n  (1.84),  and  M20  (2.05)  models.  The  red  lines  de¬ 
note  for  comparison  the  R,\  values  solely  based  on  the  con¬ 
tinuum  radiation  excluding  all  atomic  and  molecular  features. 
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Fig.  4.  Monochromatic  ta  =  1  Mira  star  radii  predicted  by  the  M  model  series  for  the  example  of  the  M16n  (model  phase  1.60),  M18n  (1.84),  and 
M20  (2.05)  models.  The  red  line  indicates  the  model  continuum  radius  excluding  atomic  and  molecular  features.  The  spectral  resolution  of  the 
model  is  0.001  gm.  Also  indicated  are  the  positions  of  H2O  and  CO  bands  after  Langon  &  Wood  (2000). 


Also  shown  are  the  positions  of  the  H2O  and  CO  bands  after 
Langon  &  Wood  (2000)  and  references  therein.  The  most  promi¬ 
nent  features  of  these  model  curves  in  the  near-infrared  region 
are  two  water  vapor  features  around  1 .4  gm  and  1 .9  /./m,  and  also 
CO  features  around  1.6//m  and  2.4pm.  The  strengths,  shapes, 
and  widths  of  these  molecular  features  depend  strongly  on  the 
stellar  phase  (and  also  on  cycle),  as  is  evident  from  the  compari¬ 
son  of  the  three  model  curves.  Also,  the  relative  strengths  of  the 
molecular  features  varies  with  stellar  phase. 

Model  M18n  provides  the  best  formal  fit  to  our  S  Ori 
AMBER  visibility  data  out  of  the  20  available  phase  and  cy¬ 
cle  combinations  of  the  M  series.  The  synthetic  visibility  values 
based  on  the  M18n  model  compared  to  our  AMBER  observation 
are  indicated  in  Fig.  1.  Here,  the  angular  photospheric  diameter 
corresponding  to  the  1.04pm  (photospheric)  model  radius  (de¬ 
fined  in  Ireland  et  al.  2004c)  is  ©MiSn  =  8.3  ±  0.2  mas,  consistent 
with  the  UD  diameter  at  near-continuum  wavelengths  1.3  pm 
and  1.7  pm  of  ©ud  =  8.1  ±  0.5  mas.  The  theoretical  Rj(ta  =  I ) 
radii  in  Fig.  4  cannot  be  compared  directly  to  the  UD  diameters 
derived  from  our  AMBER  data  (Fig.  3),  because  of  the  different 
spectral  resolution  and  because  the  model-predicted  CLVs  can 
be  very  different  from  a  UD  model  (so  that  different  radius  def¬ 
initions  are  not  equal).  The  translation  of  the  model  prediction 
into  a  UD  diameter  depends  on  the  exact  shape  of  the  bandpass- 
averaged  CLV  and  the  baselines  used.  To  compare  the  model 
predictions  to  the  measured  UD  values  in  Fig.  3,  we  fitted  UD  di¬ 
ameters  to  the  synthetic  visibility  values  of  the  M18n  model  us¬ 
ing  exactly  the  same  spectral  channels,  baseline  configuration, 
and  fit  method  as  for  our  AMBER  data.  The  resulting  model 
prediction  for  the  UD  diameter  as  a  function  of  wavelength  is 
shown  in  Fig.  3,  in  comparison  to  the  values  derived  from  the 
observation. 

Figures  1  and  3  show  that  our  AMBER  visibility  data 
can  be  described  reasonably  well  by  the  dynamic  atmosphere 
model  M18n.  The  differences  between  observations  and  model 
prediction  are  most  likely  due  to  an  imperfect  match  of  the 
phase  and  cycle  combination  between  observation  and  available 
M  models  of  the  series.  Looking  at  Figs.  3  and  4,  it  is  evident 
that  a  better  fit  to  our  AMBER  data  could  be  obtained  with  a 
model  that  shows  a  stronger  1 .4  pm  water  vapor  feature  of  the 
same  shape  compared  to  M18n  (as  seen  for  instance  in  the  case 
of  M20),  and  at  the  same  time  a  just  as  strong  but  broader  1.9  pm 
feature  (as  for  instance  in  the  case  of  M16n).  It  is  quite  pos¬ 
sible  that  such  a  combination  of  the  two  water-vapor  features 
could  appear  for  a  model  of  another  phase-cycle  combination. 
Also,  some  differences  between  M  model  predictions  and  ob¬ 
servations  of  S  Ori  are  expected  due  to  the  different  stellar  pa¬ 
rameters  of  S  Ori  compared  to  the  parent  star  of  the  M  model 
series.  Finally,  differences  can  also  be  caused  by  remaining 


uncertainties  in  the  absolute  calibration  of  visibility  values  and 
of  the  wavelength  scale.  The  S  Ori  dust  shell  of  ry  =  1.5 -2. 5 
as  modeled  in  Wittkowski  et  al.  (2007)  is  not  expected  to  have  a 
noticeable  effect  on  our  near-infrared  visibility  data,  because  its 
contribution  to  the  visibility  was  already  small  at  8  pm  and  be¬ 
cause  the  visibility  data  are  calibrated  separately  for  each  spec¬ 
tral  channel. 

In  summary,  our  AMBER  observations  of  S  Ori  generally 
confirm  the  predictions  by  the  M  model  series  and  we  find 
that  the  observed  variation  of  diameter  with  wavelength  can  be 
understood  as  the  effect  of  phase-dependent  water  vapor  and 
CO  layers  lying  above  the  photosphere.  The  M  model  series 
can  be  used  reasonably  well  to  model  the  atmosphere  of  a  Mira 
star  such  as  S  Ori  and  to  derive  a  reliable  photospheric  radius 
based  on  broadband  data.  More  such  observations  are  needed  to 
confirm  and  constrain  the  model  predictions  in  more  detail  and 
to  monitor  the  predicted  phase  dependence  of  the  strength  and 
characteristics  of  the  molecular  layers.  Simultaneously  obtained 
spectra  would  be  a  valuable  addition. 
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